LAMMPS (18 Feb 2013)
# small Peridynamic cylinder hit by projectile

units		si
boundary        s s s
atom_style      peri
atom_modify	map array
neighbor        0.0010 bin

# small target

lattice         sc 0.0005
Lattice spacing in x,y,z = 0.0005 0.0005 0.0005
region          target cylinder y 0.0 0.0 0.0050 -0.0050 0.0 units box
create_box      1 target
Created orthogonal box = (-0.005 -0.005 -0.005) to (0.005 0 0.005)
  2 by 1 by 2 MPI processor grid
create_atoms    1 region target
Created 3487 atoms

pair_style      peri/pmb
pair_coeff      * * 1.6863e22 0.0015001 0.0005 0.25
set             group all density 2200
  3487 settings made for density
set             group all volume 1.25e-10
  3487 settings made for volume
velocity        all set 0.0 0.0 0.0 sum no units box
fix             1 all nve

# spherical indenter to shatter target

variable        y0 equal 0.00155
variable        vy equal -100
variable        y equal "v_y0 + step*dt*v_vy"

fix             2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box

compute         1 all damage/atom
timestep        1.0e-7
thermo          100

#dump            1 all custom 100 dump.peri id type x y z c_1

#dump		1 all image 50 image.*.jpg type type #		axes yes 0.8 0.02 view 80 -30 adiam 0.0006
#dump_modify	1 pad 4

run             1000
Peridynamic bonds:
  total # of bonds = 335966
  bonds/atom = 96.3482
Memory usage per processor = 18.744 Mbytes
Step Temp E_pair E_mol TotEng Press Volume 
       0            0            0            0            0            0 5.0030006e-07 
     100 1.7890585e+24     552721.8            0    681881.47 1.7210968e+11 5.0030006e-07 
     200 1.3549879e+27    3097027.6            0 1.0091931e+08 1.1829272e+14 5.5130066e-07 
     300 3.3009625e+27    6331254.8            0 2.4464163e+08 1.9647265e+14 8.0862953e-07 
     400 3.815184e+27    6225081.7            0 2.8165928e+08 1.8189267e+14 1.0095118e-06 
     500 4.2580877e+27     20212686            0 3.2762196e+08 1.6249923e+14 1.2611723e-06 
     600 5.5126512e+27     30861342            0 4.2884284e+08 1.7320038e+14 1.531873e-06 
     700 1.1807414e+28     23119941            0 8.7554687e+08 2.9477434e+14 1.9278632e-06 
     800 1.2424839e+28    2407365.3            0 8.994088e+08 2.3787786e+14 2.5138992e-06 
     900 1.2358398e+28    4532419.8            0 8.9673716e+08 1.9097317e+14 3.1145903e-06 
    1000 1.2341036e+28    3219727.8            0 8.9417105e+08 1.596857e+14 3.7196039e-06 
Loop time of 5.76524 on 4 procs for 1000 steps with 3487 atoms

Pair  time (%) = 4.99867 (86.7035)
Neigh time (%) = 0.14311 (2.48229)
Comm  time (%) = 0.577355 (10.0144)
Outpt time (%) = 0.000337362 (0.00585166)
Other time (%) = 0.0457699 (0.793894)

Nlocal:    871.75 ave 910 max 835 min
Histogram: 1 0 0 0 2 0 0 0 0 1
Nghost:    1343.25 ave 1380 max 1305 min
Histogram: 1 0 0 0 0 2 0 0 0 1
Neighs:    141784 ave 170146 max 116083 min
Histogram: 1 1 0 0 0 0 0 1 0 1
FullNghs:  302269 ave 346070 max 260820 min
Histogram: 1 0 0 0 2 0 0 0 0 1

Total # of neighbors = 1209076
Ave neighs/atom = 346.738
Neighbor list builds = 46
Dangerous builds = 0
